Tuning- and order parameter in the SOC ensemble 
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The one-dimensional Oslo model is studied under self-organized criticality (SOC) conditions and 
under absorbing state (AS) conditions. While the activity signals the phase transition under AS 
conditions by a sudden increase, this is not the case under SOC conditions. The scaling parameters 
of the activity are found to be identical under SOC and AS conditions, but in SOC the activity 
lacks a pickup. 
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INTRODUCTION 

Self-organized criticality (SOC) [1 is a form of non- 
equilibrium critical behavior [2 where, in contrast to 
traditional critical phenomena, only a separation of 
timescales between external driving and internal relax- 
ation is needed for the system to develop into a scale 
invariant state. The description of SOC as the result of 
feedback between an order parameter and a tuning pa- 
rameter was put forward by Tang and Bak [2]. Later 
Dickman et al. [3 linked SOC to models with AS phase 
transitions, where a tuning parameter, the particle den- 
sity C, determines whether the system is in an active 
phase where it changes in time or in an inactive phase 
where the system is stuck in one configuration. The or- 
der parameter of these transitions is the density of sites 
about to topple, called the activity. Dickman et al. sug- 
gested that the critical behavior of SOC systems is nu- 
merically identical to that of these well understood [4 
non-equilibrium critical phenomena. 

The intuition that the two types of critical phenom- 
ena are related can be understood as follows: conserva- 
tive SOC sandpile models can be turned into AS systems 
by making the open boundaries periodic, and ceasing to 
drive the systems. Thereby, for each SOC system a corre- 
sponding AS system can be constructed. Activity in SOC 
systems transports particles in a random fashion, which 
can lead to dissipation at the open boundaries. The ac- 
tivity stops whenever the system gets stuck in an inactive 
(absorbing) state. It is then re-activated by the exter- 
nal drive. Thus, the SOC system is repeatedly driven 
through the corresponding absorbing-state phase transi- 
tion. A key assumption in the explanation of SOC in 
terms of this "AS mechanism" is that the SOC activity 
behaves just like in AS with a sudden pickup: "Slow driv- 



ing pins C at its critical value: if it exceeds £ c , activity is 
generated, and thereby dissipation" [5] • 

The main result of the present study is that this as- 
sumption is not valid. The critical properties of the AS 
transition are reflected in the finite-size scaling of the 
SOC activity. But as a function of the particle density, 
the activity shows no sign of a phase transition under 
SOC conditions. 

We note that finding the critical parameters, expo- 
nents, moment ratios etc. is incomparably easier in an 
SOC system than under AS conditions. For instance, the 
interesting long-time limit is difficult to access in simula- 
tions of AS transitions as the number of active systems in 
an ensemble shrinks exponentially with time [HI [7] ■ As a 
numerical tool for studying non-equilibrium phase transi- 
tions, SOC could be extremely useful if it was fundamen- 
tally understood (as opposed to numerically determined) 
which properties measured in SOC are universal. 

THE OSLO MODEL 

We study the Oslo model [8] in one dimension. This 
model is believed to be in the C-DP universality class 
(directed percolation coupled to a non-diffusive con- 
served field), as is the Manna model [pj. In addition 
to the results shown here, all findings were confirmed for 
the two-dimensional Oslo model and the one- and two- 
dimensional Manna models. Further details will appear 
elsewhere. 

Initialization in the AS version: G N particles for 
a fixed, given density C are added at uniformly randomly 
chosen sites in a lattice of linear size L with periodic 
boundary conditions (topology of a ring). Every site is 
assigned a threshold particle number, either 1 or 2, ran- 
domly, independently and with equal probabilities. Sites 
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exceeding this local threshold are called "active" and sub- 
ject to the relaxation process described under Toppling 
below. 

Initialization in the SOC version: The system starts out 
empty and is filled only by the External drive (below). 
Every site is assigned a random threshold as in AS. 
Toppling: Every active site is selected randomly with 
equal probability. One particle is moved to each of the 
nearest neighbors, and a new random threshold, 1 or 2, 
is assigned to the originating site. 

Avalanches: The toppling is repeated until no active site 
is left. Avalanches are defined as a sequence of topplings 
induced by a single driving step (below) . 
Microscopic time: Active sites topple at unit rate, which 
defines the microscopic time unit. In AS, observables 
are clearly time-dependent, as the system relaxes from 
the initial unstable state to complete quiescence, gener- 
ating a single avalanche. The average over initial states 
in AS is denoted (-)as- Where the time t is dropped as 
a parameter from the observables, time averaging with 
appropriate weighting has been applied. In SOC, aver- 
aged observables are denoted (-)soc- This ensemble can 
be conditioned on the density of particles £. 
Conditional activity: The activity is the density of sites 
exceeding the local particle threshold. It is the or- 
der parameter in the sense that a non-zero asymptotic 
value of lim^oo limi^oo (p a (£; L,t)) AS indicates the ac- 
tive phase. Only non-zero measurements of the instan- 
taneous activity contribute to estimates of the averages 
- in AS this is done by explicit conditioning (only ac- 
tive systems are part of the AS ensemble), while SOC 
systems are re-activated immediately whenever they fall 
into an absorbing state. Nonetheless, it is important to 
distinguish between 

1) the AS activity (p q (C;^))as measured at constant 
(conserved) externally set ( and 

2) the SOC activity (p a (& ^))soc conditioned on the 
value £ of the fluctuating density, £(t), in the SOC en- 
semble. 

Boundary conditions: In the AS version, boundaries are 
periodically closed and no dissipation of particles takes 
place. In the SOC version, boundaries are open, i.e. top- 
pling across a boundary removes particles from the sys- 
tem. 

External drive: The AS models are not driven. Time se- 
ries (see below) terminate when all activity has ceased. 
This happens, in a finite system, with a rate that is 
bounded away from zero, i.e. in any finite system the 
avalanche eventually comes to a halt (below a trivial den- 
sity limit) . While the model in the AS version is restarted 
by resetting it to a new initial condition, the SOC version 
is driven externally to compensate the loss of particles at 
the boundary. The driving occurs only at quiescence, 
thereby completely separating the timescales of toppling 
and driving. 

Time series in the AS version: An ensemble of systems 



is prepared, and observables are averaged over the mem- 
bers of the ensemble (of systems that are still active) at 
equal times, producing time series of fluctuating observ- 
ables. We focus here on moments of the activity in the 
quasi-stationary state, (Pa(C)) AS ; that is, moments of 
the activity after a transient but conditional to activity 

DIED- 

Time series in the SOC version: Because the SOC sys- 
tems are driven externally, they do not need to be reset 
to the initial state and restarted. After a sufficient tran- 
sient, which is dominated by the time it takes the exter- 
nal drive to fill the system to a density near £f OC > t ne 
observables fluctuate about an average value. 

Typical SOC observables include the avalanche size, 
defined as the number of topplings, and the duration, de- 
fined as the number of microscopic time units that elapse 
during an avalanche. 

Studies of AS systems typically focus on stationary and 
dynamic properties more commonly studied in critical 
phenomena, such as the finite-size scaling of the order 
parameter and its fluctuations, survival probabilities, and 
spreading of activity. These observables will be referred 
to as AS observables. 



SIMILARITIES BETWEEN AS AND SOC 

Given the claim that corresponding observables scale 
identically under SOC and AS conditions, surprisingly 
few direct comparisons exist. Specifically, we are not 
aware of any studies of critical scaling (in the tuning pa- 
rameter) of observables under SOC conditions. However, 
it has been reported that 

• the asymptotic (L — > oo) particle density in an 
SOC system coincides with the critical density of 
the corresponding AS system [TT1 [T2l [H] , 

• the avalanche sizes measured in an AS system, that 
is, the response of a finite quiescent AS system at 
the critical particle density to the addition of a par- 
ticle, obey the same finite-size scaling as avalanche 
sizes in bulk-driven SOC systems [TTJ [TJ] . 

• scaling relations between avalanche exponents de- 
scribing SOC observables and exponents character- 
ising AS systems [TTJ [T31 [T3J [T5] are apparently 
valid. This implies that (p a (L)) soc follows the 
same scaling as (p a (Cc, L)) AS . 

DIFFERENCES BETWEEN AS AND SOC 

These similarities are remarkable because the ensem- 
ble sampled by SOC systems is very different from that 
sampled by AS systems. 
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Our main finding is illustrated in fig. [T] Appropri- 
ate rescaling with L of (p q (C;^))as anc ^ C collapses the 
curves (p a (Ci -^))as- The same scaling parameters pro- 
duce a collapse of the curves (p a (Ci -^))soc ^ e * w0 
collapsed curves differ. The AS activity hence follows 



(Pa(C,L)) AS = 



.AS 



L 



-p/u ±n AS fn AS 



G AS (a2 a ({ - QL 1 '^) (f) 



where a A ^ denotes metric factors [T7] , G AS a scaling func- 
tion and v± and /3 are the usual [4] critical exponents. 
The SOC activity follows 

<Pa(0 L)) soc = af Oi-^^floo^so c (c _ 

(2) 

with numerically identical constants as in eq. (|2]). Cru- 
cially, however, the scaling functions Q soc and G AS are 
fundamentally different, with (p (C;^))soc having no 
powerlaw pickup. In other words, the ensemble probed in 
the SOC regime lacks the one key feature that motivates 
the AS mechanism. 

Since the scaling parameters are identical under SOC 
and AS conditions, our findings do not contradict the 
many observations of valid scaling relations between SOC 
and AS exponents. But our findings call into ques- 
tion the physically appealing narrative of the AS mech- 
anism, where an SOC system is pushed through an AS 
phase transition, into the active phase by a slow drive 
and strongly repelled from there as a sudden increase 
in activity above £ c leads to dissipation. This sud- 
den increase does not exist in SOC. Our critique [18] 
was based on the disbelief that a linear feedback loop 
like the equation of motion eq. ^ could drive a non- 
linear system like the Oslo model to a (highly non-linear) 
phase transition. Fig. [I] specifies the discord between 
SOC and AS ensemble, invalidating the key proposition 
(i°o)soc (C; L ) = (pa) as (C; L ) underlying the AS mecha- 
nism. 

Neither the order parameter nor its variance (not 
shown) have special features near the critical particle 
density Q c . Whereas a non- linear ^-dependence is clearly 
seen under AS conditions, (p a (C'i L)) AS , over the range of 
( sampled by the SOC system, (p a (C, L)) soc increases 
approximately linearly with £. 

The distribution of the SOC system's tuning parame- 
ter P(Csoc(*); L) is well described by a Gaussian. Fig. [2] 
shows its width tTg OC (£) , which is a measure for the ap- 
parent critical region in the SOC model. It is found to 
shrink as L^ x ^ u± , which is the minimum rate required for 
universal finite size scaling in the light of the AS mecha- 
nism |18j . 



DISCUSSION 



The activity (pa(O)soc m tne SOC model does not 
signal the onset of a phase transition around £ c by a kink. 




FIG. 1: Rescaled order parameter versus rescaled reduced 
density (circles: L = 256, squares: L — 1024, triangles: L = 
4096). The dependence on £ of the SOC order parameter 
{p a } soc (Ci L) (open symbols), looks very different from the 
AS order parameter {p a ) AS (Ci L) (filled symbols). The same 
scaling parameters, f3 = 0.26, u± — 1.33 [19], and f c = 1.73260 
[14| are used for SOC systems and for AS systems. The active 
phase can be seen in AS mode as a power-law pick-up in the 
scaling function but is invisible in SOC mode. For the extreme 
values of £ reached by the SOC system, some error bars are 
curiously small, which is due to small sample sizes at extreme 
parameter values. 
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FIG. 2: The width of the distribution P((soc(t); L) as a func- 
tion of system size L shrinks as L^ 1 ^ 1 - , enabling the observa- 
tion of universal (AS-) exponents in the SOC ensemble. The 
value v± — 1.33 was used for the rescaling in this figure. The 
re-scaled quantity varies by less than 4% between L — 256 
and L = 65, 536. 



This seems to contradict the widely accepted explanation 
of how the SOC system finds the critical point. 

The comparatively suppressed activity for £ > £ c under 
SOC conditions confirms that the difference between the 
ensembles is significant: In order for the SOC system to 
advance deep into the high-£ region, low dissipation and, 
presumably, low activity are needed. The particle density 
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£ under SOC conditions thus has a different effect from 
£ under AS conditions. 

Scaling laws relating SOC and AS exponents are 
based on the assumption that moments (/9^(L)^ soc and 
(p^(C c ;£)) AS follow the same scaling. The validity of 
such scaling laws has been confirmed repeatedly in the 
literature [T3J [2U] . We have added to these studies the di- 
rect observation of identical scaling of (p a (C; L)) soc and 
(p a (£; L)) AS , i.e. over a range of £, but we have also 
shown where the equivalence between corresponding AS 
and SOC models breaks down. 

Given the clear difference between the AS and SOC- 
behavior of such a defining observable as the order pa- 
rameter near criticality - why do we observe identical 
finite-size scaling and, even more puzzlingly, why iden- 
tical critical densities? The ideal theory answering this 
question includes a recipe for turning any phase transi- 
tion into SOC - if that is only possible for some "special" 
AS transitions (maybe there is only one such universality 
class), the theory has to explain this restriction, and clar- 
ify how much of our understanding of phase transitions 
applies to AS transitions at all. 

The AS mechanism of self-organization, as formulated 
by Dickman et al. can be paraphrased in a mean-field 
equation of motion for the tuning parameter [3] , 
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d t at) = h(t,L)-e(t,L)p b a ^(t,C,L) 



(3) 



where h(t; L) is the rate of increase in £, measured on 
the microscopic time scale, due to the addition of par- 
ticles, and e(t; L)pa° c (t,(', L) is the loss rate. Even if 
one accepts that such a simple feedback loop can drive 
a system into a phase transition, the models we inves- 
tigated (Oslo and Manna 1-d and 2-d) do not seem to 
be governed by its supposed equation of motion, eq. Q3J) , 
because the detailed behavior of p^° c in this equation 
cannot be equated to that of p^ s . 

Because it is so widely accepted in the literature 
that SOC probes a corresponding AS phase transition 
[5J [TTJ [TSJ, |5T] , it is often difficult to see which results are 
based on this assumption and which support it indepen- 
dently. We emphasize that the scaling relations between 
AS and SOC are not an explanation of the equivalence 
between AS and SOC; their apparent validity constitutes 
the observation of this equivalence that still needs to be 
understood. 
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